Conductance peaks in open quantum dots 
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We present a simple measure of the conductance fluctuations in open ballistic chaotic quantum 
dots, extending the number of maxima method originally proposed for the statistical analysis of 
compound nuclear reactions. The average number of extreme points (maxima and minima) in the 
dimensionless conductance, T, as a function of an arbitrary external parameter Z, is directly related 
to the autocorrelation function of T{Z). The parameter Z can be associated to an applied gate 
voltage causing shape deformation in quantum dot, an external magnetic field, the Fermi energy, 
etc.. The average density of maxima is found to be (pz) = az/Zc, where az is a universal constant 
and Zc is the conductance autocorrelation length, which is system specific. The analysis of (pz) 
does not require large statistic samples, providing a quite amenable way to access information about 
parametric correlations, such as Zc- 
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Introduction - The statistical properties of the elec- 
tronic transport in ballistic open quantum dots (QDs) 
have been intensively studied over the last decades ■ 
In such systems, the conductance can be described by 
the Landauer formula and, for QDs containing a large 
number of electrons, the random matrix theory (RMT) 
provides an excellent statistical description of the under- 
lying chaotic electronic dynamics at the Fermi energy 
[3|, |4| . RMT explains the observed universal conductance 
fluctuations in QDs, which depend only on the QD sym- 
metries, such as time-reversal, and on the number of open 
modes TV connecting the QD to its source and drain reser- 
voirs [ij. 

In the semiclassical limit of large N, the transmission 
statistical fluctuations are accurately modeled by Gaus- 
sian processes. In practice, it has been experimentally 
observed Q and theoretically explained [1] that, even for 
small values of N and at very low temperatures, dephas- 
ing quickly brings the QD conductance fluctuations close 
to the Gaussian limit. 

The conductance in open ballistic QDs exhibits ran- 
dom fluctuations as an external parameter, such as a 
magnetic field B or an applied gate voltage Vg, is var- 
ied. By identifying running averages with ensemble aver- 
ages, it is customary to accumulate statistics by varying 
as many parameters as the experimental set up allows. 
This invites one to ask whether useful statistical informa- 
tion can be extracted from the analysis of a single con- 
ductance curve. Inspired by the formal analogy between 
conductance and compound-nucleus Ericson fluctuations 
Q we show that the answer is positive. More specifi- 
cally, we further develop ideas originated in the context 
of nuclear physics [S], to calculate the conductance av- 
erage density of maxima and show its relation with the 
conductance autocorrelation function. As a result, we 
propose a new universal measure for the conductance of 
ballistic open QDs. 



Theoretical framework. - We consider the standard 
setting of a two-probe open quantum dot coupled by leads 
to a source and a drain electronic reservoirs. We also 
assume that the source (drain) reservoir is coupled to the 
quantum dot by a lead that has iVi (TV2) open modes. 
The scattering matrix S describing the electron fiow is 
given by Q 

s^{y,) (1) 

where r (r') is the A^i x iVi {N2 x N2) matrix containing 
the refiection amplitudes of scattering processes involving 
channels at the source (drain) coupled leads, while t {t') 
is the iVi X A'2 {N2 x Ni) matrix built by the transmis- 
sion amplitudes connecting channels that belong to the 
source-coupled lead to the modes at the drain-coupled 
lead (and vice- versa). 

At zero temperature, the linear conductance G of an 
open quantum dot is given by the Landauer formula 

G = T with T = tr(th) (2) 

h 

where the factor 2 accounts for spin degeneracy and T 
is the dimensionless conductance or transmission, which 
typically depends on A^i and N2, the quantum dot shape, 
the external magnetic field B, the electron energy e, etc.. 

In the limit of large number of open modes, the average 
transmission for a chaotic QD is Q 

(T) = — (3) 

^ ' N1+N2 4 ' ^ ' 

where (• • •) indicates that an ensemble average was taken 
and 13 — 1 — 2) corresponds to the orthogonal (uni- 
tary) case of preserved (broken) time-reversal symmetry. 
In the same limit, the transmission correlation function 
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reads M-h 



(T«(e,X)rft(e',X')) 



var^(T) 



[l + (5X/X,)2)]2 + (fc/r)2' 



(4) 

where var/3(T) = (1 + (5^^i)/16. To simplify the nota- 
tion, we introduce T*^{e,X) = T{e,X) - (T), where e is 
the electron energy and X is a generic parameter that 
describes a certain quantum dot shape belonging to a 
path of deformations caused by, for instance, applying 
a certain gate potential. The correlation function given 
by Eq. 2] is universal, with correlation length scales, Xc 
and r, that are system dependent. There is a simple 
expression that relates F to the mean resonance spacing 
A, namely, T = (iVi + 7V2)A/27r The correlation 

length r is generally different from the "lifetime" or de- 
cay width, which is twice the imaginary part of the pole 
energy of the scattering matrix. Both quantities only co- 
incide when r ^ A, a condition never met in open QDs. 

Density of maxima. - The transmission or dimension- 
less conductance T{Z) as a function of a generic param- 
eter Z (either s or X) has a maxima in the interval 
[Z, Z + SZ] if 

T'{Z) > and T'{Z + 6Z) <0 , (5) 

provided 6Z is small. In this case, Eq. [5] implies that 



T"{Z)5Z > T'{Z) > 



(6) 



For convenience we introduce T' and T" to denote the 
first and second derivatives of the dimensionless conduc- 
tance T with respect to Z. 

The joint probability distribution P{T', T") allows one 
to obtain the average density of maxima {pz) The 
probability to find a maximum in the interval [Z, Z -\- 6Z] 
is the integral of P over the region defined in Eq. |6l that 
is 



dT" / dT'P{T',T 

'OO Jo 



I TJlrril rTill\ 



-5Z 



dT"T"P{Q,T") EE 6Z {pz) 



(7) 



Let us infer P{T',T") by examining the lowest mo- 
ments of T' and T". Since the statistical properties of the 
dimensionless conductance are invariant under Z trans- 
lations, T' and T" have zero mean. Their variance is 
directly related to the correlation function 



Czi6Z) = {T^{Z + a6Z) T^{Z - bSZ)) 



(8) 



which does not depend on the choice of a and b, provided 
a+b = 1. Neither do so the derivatives of with respect 
to SZ, which leads to 



{[Tf) 
{TT") 

{[T'T) 



SZ=Q 



SZ=0 



sz=o 



(9) 



and {TT') = {T'T") = 0. These results coincide and 
expand those obtained in Ref. [13. 

We use the above relations and the maximum informa- 
tion principle to built the joint probability distribution 
of the transmission T and its derivatives, T' and T" . The 
distribution P{T', T") is found by integrating over T, and 
gives 



p(o,r") = 



1 [T"]2 
expl --■ 



V{[T'?) {[T'T) 2([r"]2 

Thus, the integral in Eq. [7] renders 



(10) 



^^^^ 2^V ([T']2) 



(11) 



This result, obtained with the help of the maximum infor- 
mation principle, is expected to be accurate in the large 
iVi -I- N2 limit due to the central limit theorem [l^ . In the 
opposite limit of small N1 + N2, the ratio (T) /[var(r)]i/2 
is no longer large and the constraint T > gives raise to 
non trivial correlations between the transmission and its 
derivatives, as well, as deviations from the Gaussian dis- 
tribution. 

In the case where the external parameter Z is the elec- 
tron energy e and Ni+ N2 ^ 1, the correlation function 
given by Eq. [5] reduces to a Lorentzian 



varp (T) 



Such correlation function gives 



0.55 



(12) 



(13) 



Hence, by counting the average number of maxima in 
conductance one can infer the conductance correlation 
width. This idea was originally proposed as comple- 
mentary to the analysis of the Ericson fluctuations in 
compound nucleus reactions [l^ [161. The analysis of 
Ref.hjl seemingly different from ours [17| , gives the same 
result as above. 

Support to our analytical findings is provided by nu- 
merical simulations employing the Hamiltonian approach 
to the statistical iS-matrix jl8j , namely 



S{e) 1 - 2TTiW^e -H + iirWW^'y^W 



(14) 



where e is the electron propagation energy and H is the 
matrix of dimension M x M that describes the reso- 
nant states. H is taken as a member of the Gaussian 
orthogonal (unitary) ensemble for the (broken) time- 
reversal symmetric case. The matrix W of dimension 
M X (A^i + N2) contains the channel- resonance coupling 
matrix elements. Since the H matrix is statistically in- 
variant under orthogonal (/3 = 1) or unitary (/? = 2) 
transformations, the statistical properties of S depend 
only on the mean resonance spacing A, determined by H, 
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and W'^W . We assume a perfect coupling between chan- 
nels and resonances, which corresponds to maximizing 
the average transmission following a procedure described 
in Ref. [18| . In this paper we restrict our numerical anal- 
ysis to the /? = 1 case and, for simplicity, we take the 
case of = A^i = N2. We benchmarkcd the accuracy 
of the simulations by an extensive comparison between 
numerical simulations and analytical results [9] for (T) 
and var(T) as a function of N. 

Figure [1] illustrates the transmission T{e) for a typical 
realization of the matrix model given by Eq.ll4| for N — 5 
perfectly coupled modes close to the band center. The 
transmission correlation length is given by the Weisskopf 
estimate namely, T = {Ni + N2)A/2n = iVA/7r. 
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FIG. 1: Typical dimensionless conductance T as a function 
of e for N — 5 and perfectly transmitting modes (no direct 
processes). Black dots stand for the numerical results for a 
single realization of H and the dotted line indicate the RMT 
prediction for (T). 

Figure [5] shows the transmission autocorrelation func- 
tion Ce((5e) obtained from the model given by Eq. [141 
for N perfectly coupled modes. The ensemble average 
is taken over 200 realizations of the H matrices with 
M = 200. The random matrix theory [l^ predicts an 
autocorrelation length F = NA/tt, which is nicely veri- 
fied by the simulations. 

Figure [3] shows the average density of maxima (pe) in 
units of 1/r as a function of the number of open chan- 
nels N. We observe that the agreement with the Gaus- 
sian process prediction becomes remarkably good as TV 
is increased. 

Let us switch our analysis to the case where an external 
parameter modifies the QD Hamiltonian, namely, H — 
H{X). Taking Se = 0, the transmission autocorrelation 
function, Eq. |4l becomes a Lorentzian squared [l0| 



CxiSX) 



var^(T) 



[1 + {6x/x,y 



(15) 
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FIG. 2: Normalized transmission autocorrelation function 
Ce{Se) = Ce(5e)/var(r) as a function of the energy 5e. Sym- 
bols correspond to ensemble averages for different number of 
channels A^. The statistical error are smaller than the symbol 
sizes. The solid line stands for Cs{Se) given by Eg. 1121 




FIG. 3: Density of maxima {pe)r as a function of the number 
of open channels A'^. The symbols with statistical error bars 
correspond to our numerical simulations. The dashed line 
stands for the Gaussian process prediction. 



This correlation function gives for (px) 

, , 3 0.68 

{px) 



7rV2Xc 



Xr. 



(16) 



The above result is new and is tested through numerical 
simulations in what follows. 

To statistically model H{X) we take H = Hi cos X + 
H2 sin X , where both Hi and H2 belong to a Gaus- 
sian ensemble. The transmission T{X) is obtained by 
computing the S matrix defined by Eq. [T3] at e = for 
1,000 realizations of the H{X) with M = 200. Figured 
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FIG. 4: Normalized transmission autocorrelation function 
Cx{&X) = Cx{5X) /Y&iiT) as a function of the parameter 
8X/Xc- Symbols correspond to numerical simulations for dif- 
ferent A'^. The statistical error bars are smaller than symbol 
sizes. Solid line is given by theory, Eq. 1151 Insert: Xc versus 
N^^'^ showing a linear behavior, as indicated by the solid line. 



shows that a Lorentzian squared adjusts very nicely the 
numerically obtained correlation functions upon rescal- 
ing X hy Xc. 

In distinction to the previous case, where a simple an- 
alytical expression for T is known, here we determine Xc 
numerically. Using semiclassical arguments, it can be 
shown [21] that the effect of a perturbation grows dif- 
fusively with the electron dwell time in the quantum 
dot, which scales as ^ 1/A^. Hence Xc ^ ^fN , in ex- 
cellent agreement with our numerical findings, shown in 
the inset of Fig. S) 

Figure [5] summarizes our numerical results for the case 
of parametric Hamiltonian changes: The density {px) in- 
creases with N and rapidly saturates at a value in very 
good agreement with our Gaussian process prediction 
given by Eq. (I16p . Conclusions- In this work we extended 
the number of maxima method, originally employed in 
compound nuclear reactions, to open chaotic QDs. We 
have shown that the average density of maxima in the di- 
mensionless conductance is inversely proportional to its 
autocorrelation length. For parametric variations that 
give rise to a Lorentzian-like transmission autocorrelation 
function, like variations in the electron energy, the uni- 



versal proportionality constant is y/Z/ir. For parametric 
changes that lead to squared Lorentzian-like transmission 
correlations, such as gate potential variations, the univer- 
sal proportionality factor is 3/(7r-\/2). These results are 
obtained by assuming that the transmission derivatives 
are Gaussian distributed, which is expected to be rather 
accurate in the semiclassical limit of large N . We em- 
ploy numerical simulations to infer the precision of our 
results for an arbitrary N . We show that even for moder- 
ate values of N the semiclassical prediction gives already 




FIG. 5: Density of maxima (px) as a function of the number 
of open channels A'^ in units of Xc. The symbols with statis- 
tical error bars correspond to our numerical simulations. The 
dahed line stands for the theoretical prediction, Eq. [T6] 



qualitatively good results, within about 10% precision. 
Our results may prove useful for the analysis of measure- 
ments of the transmission in chaotic quantum dots: By 
counting the maxima of a simple magnetoconductance 
trace, it is possible to estimate with a rather good pre- 
cision the dimensionless autocorrelation function. More 
generally, ballistic mesoscopic systems (and potentially 
diffusive ones) showing conductance fluctuations, such as 
graphene flakes |22h24| . are also potentially amenable to 
this analysis. 
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